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We investigate the effect of impurities in multiferroic materials using an equation of motion 
approach for the spin dynamics of the host multiferroic compound. We model the impurities as 
a two-level system and focus on the regime where the impurity spins relax slowly. When the 
impurity strength is weak the host spins oscillate with no decay and the electric polarization is not 
affected. However as the impurity strength is increased the host spin components get damped and 
the electrion polarization is suppressed. Since polarization in multiferroic materials is driven by 
magnetic ordering we conclude that the presence of impurities is detrimental to multiferroicity. 



PACS numbers: 75.80.+q, 75.47.Lx, 77.80.-c 



Ferroelectricity in multiferroics is driven by magnetic 
ordering. The non-collinear spin arrangement below the 
transition temperature drives the system to a space- 
inversion and time-reversal broken symmetry state [H, 0] • 
This causes an electric polarization to be generated. The 
electric polarization in these materials can be controlled 
by a magnetic field and the magnetization by an electric 
field, but are there any factors which are detrimental to 
this controllability? In this Brief Report, we answer such 
a question by studying the effect of impurities in multi- 
ferroic compounds. Present theoretical [3, 0, H 0, Q 
and computational treatments || of multiferroic systems 
have not considered the effect of impurities. 

A multiferroic system with impurities can be described 
by the microscopic Hamiltonian, H , which generates the 
spin-spiral ground state configuration and includes a cou- 
pling of the host spins with the impurity spins 



H = H, 
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where H m f is the multiferroic Hamiltonian [ij and in- 
cludes the following terms 



H m f = H ex + Hf-in + H 
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In the above Hamiltonian H ex is the Heisenberg exchange 
Hamiltonian for nearest-neighbor spin interaction, Si and 
Si+i, of the host multiferroic material. The exchange 
parameter J between nearest-neighbor spins is consid- 
ered to be constant. The term Hum represents the 
Dzyaloshinskii-Moriya interaction and the spin-lattice in- 
teraction strength is denoted by A. The dynamics of the 



displacement at site i and the canonically conjugate mo- 
mentum is P{. In Hkin, K and M are the spring constant 
and the effective mass of Ui respectively. The easy plane 
spin-anisotropy term is given by the last expression, H an , 
where the parameter D > 0. The easy plane in this model 
is the x — z plane. 

We model the impurities as a two-level system de- 
scribed in terms of an effective spin- 1/2 The impu- 
rities are assumed to be distributed uniformly through- 
out the host material and couple with the spins of the 
multiferroic compound via an exchange interaction. It 
is assumed that the impurities are diluted enough that 
they do not form their own fully magnetized sublattice 
system. We therefore neglect the spin-spin interaction 
between the impurities. The Hamiltonian, Hf®, of the 
impurity-multiferroic exchange interaction between the 
impurity spins and the neighboring host atoms is 
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where the impurity spin components at R m is denoted 
by s(Rm) = (s Xm (Rm),s Vm (Rm),s Zm (Rm)) and the spins 
of the multiferroic system located near the m-th impu- 
rity is given by S(R m + r v ) = (S x (R m + r v ), S y (R m + 
r u ),S z (R m + r v )). The local co-ordinates for the fri- 



th impurity spin are a m = x m , y n 



and a = x,y,z 



are the co-ordinates of the host spin components. The 
anisotropic exchange parameters of the interaction be- 
tween s(R m ) and S(R m + r u ) are A am a (R m ,f v ) where 
r v is the position vector from the m-th impurity spin to 
the i-the spin, S^. 

To obtain a quantum mechanical version of Hf^ p we 
express the impurity spin components in terms of the 
creation (cTj and the annihilation operator (c m ) respec- 
tively. Using canonical commutation relations in Eq. [7] 
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we can derive the impurity Hamiltonian, Hi mp , as 

Hi m p = E &m(t = 0)n m 4- (8) 
HZ P = E ^ • fr™ K - 1/2) + e. m , ym (4) 

z,m 

+ JM] (9) 



i/ m / part of the Hamiltonian and find 

h\ = 4 JSf cos qa 

■sin 2 qa(Sf +1 + Sf_i) 
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and e«. 



where f2 m (t = 0) represents the initial splitting be- 
tween the energy levels of the two-level system used to 
model the impurity spins, n m is the population of the 
impurity spins at that site for a given temperature, T. 
The coupling coefficients are given by the expressions 

^Y,u A ^ m ,a(Rm 1 r„). The coefficient 
£x m ,a;y m ,a is related to the fast relaxation process and 
£z m ,a is for the slow relaxation process. 

The fast relaxation mechanism accomplishes direct ex- 
citation of the impurity ions by the magnetic oscillations 
of the system. The slow relaxation mechanism is based 
on the modulation of the impurity energy levels by the 
oscillations of the magnetic system and is the relaxation 
process considered in this paper [To| . At present there is 
no experimental study in multiferroic compounds which 
investigates the effects of impurity on spin relaxation rate 
and its consequences on electric polarization. The con- 
sideration of slow relaxation as a mechanism for spin re- 
laxation in multiferroics is therefore an assumption in 
this paper. 

We use Eqs. (2) - (6) and (7) - (9) to study the effect 
of impurities on the spin relaxation. For this purpose we 
implement an equation of motion approach for the spin 
dynamics [l(| El of the system. Consider the equation 



dt 
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where 7 is the gyromagnetic ratio and the effective field, 
h e ff, is given by 



h eff ~ 



dH 

dS t 
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The second term in the above expression arises from the 
impurity spins. To obtain the effective field contribution 
we introduce a rotating local co-ordinate system (£, 77, () 
and the effective field hi = (h^h^hj) in this frame of 
reference 1121. We first obtain the contribution from the 



where q is the wavector and a is the lattice spacing. 

The contribution to the effective field from the impu- 
rity part of the Hamiltonian can be obtained by consid- 
ering the kinetics of the m-th impurity spin population 
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where rii m (n m ) is a relaxation coefficient for the im- 
purity spins. At each instant of time the population 
n m = ( c ln c m) relaxes to the equilibrium value ny(f2 m ) = 
[exp(ri m /fcsT)-|-l] _1 corresponding to the dynamic split- 
ting fi m = Q m (t = 0) + SQl n (t) where 5Q^(t) = e Zm ■ S t . 
For the case of slow relaxation the impurity spins mod- 
ulate slowly, I^OjyJ/fcgT << 1, and the impurity contri- 
bution to the effective magnetic field is Q 
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The third term in the above equation, Eq. 1161 is responsi- 
ble for generating relaxation and the first and the second 
term introduces an anistropy. We will neglect the effect 
of the anisotropies generated by these terms and focus 
on the effect of the relaxation term. 

Using Eq. [16] and the equation of motion (Eq. [TOf con- 
tribution for the spin components in the rotating frame, 
(Sf,S^,S^), the effective magnetic field can be written 
as 



dS[ 
dr 

dSl 
dr 

dS C 
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where r 



(11) ^£ = 2 ( D + 2 J cos 2 SfS? - >yf3sA 



dT 

fit. The damping co-efficient (3 

cxp(n /fc B T) 
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r l|feB T; [cxp(n /fc B T)+i]^ where k B * the Boltzmann 
constant. The impurity concentration, c imp = "fn imp z imp 
where n imp is the number of impurities per unit volume 
and Zi mp is the average number of magnetic neighbors for 
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there is no polarization along the 6-axis, so the expecta- 
tion value of the displacement field along that direction 
is zero, (u y ) — 0. We therefore ignore the contribution of 
these terms in the derivation of the equation of motion. 
In the equation of motions, Eqs. [T7l-fT9l we have ignored 
terms greater than second order in the spin components. 
Also, the effective magnetic field terms which are pro- 
portional to the A coupling do not appear in the final 
expressions as these terms belong to a fourth order spin 
contribution and are neglected. Therefore, in the present 
theory (second order in the spin components) the spin- 
lattice interaction has no role to play. 

We solve Eas.[T7l-fT9lfor appropriate initial conditions 
to obtain Fig[T] To begin with we consider the spiral in 
the x — z plane. Parameters relevant for the multiferroic 
material TbMn03 [13] are used for the simulation: J = 
0.15 meV, and D — 0.4 meV. We choose qa — ir/4 and 
7 = 2. From Fig. [1] we observe that as the contribution 
of the /3— parameter is increased their is a pronounced 
damping of the S? and spin components (see Fig. [IJ . 
In the presence of damping the S^ component reaches an 
equilibrium value and its oscillations are suppressed also. 

Electric polarization, P, in multiferroic materials can 
be generated by a spin current, j s a Si x Sj [4]. The 
polarization vector defined as P oc eV, x (Si X Sj), implies 
that the spin current for the spin arrangement chosen 
in this calculation should be along the y— axis and the 
polarization along the a;— axis. Using the expression for 
spin components in the rotating frame [12] we have for 
the magnitude of the x— component of P 



(c)/3 = 10 



\P x \ = [S 2 + (S^) 2 } S m(A6) 
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FIG. 1: (Color online) Relaxation of spin components, 
(S i ,S t >,S c ), for various values of the damping coefficient (3. 
The range for /3 is varied from 0.1 to 10 by a factor of ten. 
The spin component, 5% responsible for generating the elec- 
tric polarization (see Eq. 120 along the a;— axis is damped 
leading to a suppression of the induced ferroelectricity. 



one impurity. The damping coefficient is dependent on 
the impuritiy concentration (ci mp ), the temperature (T), 
the energy scale of the impurities (fl m (t — 0)) taken to 
be a constant value f2 m (t = 0) = Q , and the impurity 
coupling coefficient e m — (e, 0,0). 

In deriving the spin dynamics of the multiferroic sys- 
tem we have considered the following facts. The original 
effective magnetic field generates terms which couple the 
u y displacement field. However, in the polarized state 



where we have used the fact that Sf +1 rj Sf and set 
either component to S^, A6> = 9i — ftj+i is the angular 
difference between the spins in the nearest neighbor sites. 
In the above equation we have replaced the S^ compo- 
nent of the spin with an equilibrium value, S, as seen in 
Fig. [U From Eq. [3U] we observe that as the damping co- 
efficient is increased, a suppression of the S'' component 
(see Fig. [T]) occurs leading to a suppression of the electric 
polarization, \P X \. 

Based on the above analysis we conclude that impu- 
rities are detrimental for multiferroic materials. Impuri- 
ties lead to a damping of the spin polarization compent 
amplitudes and in turn to a suppression of ferroelectric 
polarization. 

TD acknowledges helpful discussions with Kingshuk 
Majumdar. 
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